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Abstract 

The relative production rate of B® and B° mesons is determined with the hadronic 
decays B® — > D~ir + and B° — > D~ K + . The measurement uses data corresponding 
to 1.0 fb _1 of pp collisions at a centre-of-mass energy of ^/~s = 7 TeV recorded in the 
forward region with the LHCb experiment. The ratio of production rates, f s /fd, is 
measured to be 0.238 ±0.004 ±0.015 ±0.021, where the first uncertainty is statistical, 
the second systematic, and the third theoretical. This is combined with a previous 
LHCb measurement to obtain f s /fd = 0.256±0.020. The dependence of f s /fd on the 
transverse momentum and pseudorapidity of the B meson is determined using the de- 
cays B® — > D~ir + and B° — > D~ir + . There is evidence for a decrease with increasing 
transverse momentum, whereas the ratio remains constant as a function of pseudora- 
pidity. In addition, the ratio of branching fractions of the decays B° — > D~ K + and 
B° -> D-tt + is measured to be 0.0822 ± 0.0011 (stat) ± 0.0025 (syst). 
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1 Introduction 



The ratio of fragmentation fractions f s /fd quantifies the relative production rate of 
B° s mesons with respect to B° mesons. Knowledge of this quantity is essential when 
determining any B® branching fraction at the LHC. The measurement of the branching 
fraction of the rare decay B® — > [l] is the prime example where a precise measurement 

of f s /fd is crucial for reaching the highest sensitivity in the search for physics beyond the 
Standard Model. The branching fractions of a large number of B° and B + decays have 
been measured to high precision at the B factories (2] , but no B® branching fraction is yet 
known with sufficiently high precision to be used as a normalisation channel. 

The relative production rates of b hadrons are determined by the fragmentation fractions 
fui fdi fs, fc aud /a, which describe the probability that a b quark will hadronize into a 
B q meson (where q = u, d, s, c), or a b baryon, respectively]^] The ratio of fragmentation 
fractions f s /fd has been previously measured at LHCb with hadronic (3) and semileptonic 
decays (3] , and the resulting values were combined [i] . 

In this paper, the ratio of fragmentation fractions f s j fa is determined using B° s — > D~tt + 



and B — > D~K + decays collected in pp collisions at a centre-of-mass energy of y/s = 7 TeV, 
with data corresponding to an integrated luminosity of 1.0 fb~ recorded with the LHCb 
detector. Since the ratio of branching fractions of the two decay channels is theoretically 
well understood [5], their relative decay rates can be used to determine the ratio of 
fragmentation fractions for B° s and B° mesons through 

f, B(B°^D-K + )e DK N Ds7r 
f d B(B^D-n+)e Dsn N DK 

= $PS 
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(1) 

where N corresponds to a signal yield, e corresponds to a total efficiency, t b o/t b o = 
0.984 ± 0.011 (6] corresponds to the ratio of lifetimes and B(D~ -> K + n~n-) = (9.14 ± 
0.20)% |7] and B{D~ -> K + K~n~) = (5.50 ± 0.27)% |8| correspond to the £>^ meson 
branching fractions. The factor M a = 1.00 ± 0.02 accounts for the ratio of non-factorizable 
corrections |9|, Mf = 1-092 ± 0.093 for the ratio of B? g s — > form factors [10 , and 
$ps = 0.971 for the difference in phase space due to the mass differences of the initial 
and final state particles. The numerical values used for the CKM matrix elements are 
\V US \ = 0.2252, \V ud \ = 0.97425, and for the decay constants are /„ = 130.41 MeV, 
fx = 156.1 MeV, with negligible uncertainties, below 1% 12]. The measurement is not 
statistically limited by the size of the B° — > D~K + sample , and therefore the theoretically 
less clean B° — > D~tt + decays, where exchange diagrams contribute to the total amplitude, 
do not contribute to the knowledge of f s /fd ■ 

The ratio of fragmentation fractions can depend on the centre-of-mass energy, as well 
as on the kinematics of the B?, meson, as was studied previously at LHCb with partially 
reconstructed B decays 12]. The dependence of the ratio of fragmentation fractions on 



1 Charge conjugation is implied throughout this paper. 



1 



the transverse momentum and pseudorapidity 77 of the B9s meson is determined using 
fully reconstructed B° — > D~n + and B® — > D~n + decays. Since it is only the dependence 
that is of interest here, the more abundant B° — > D~n + decay is used rather than the 
B° — > D~K + decay. The B° — > D~K + and 5°-> D~tt + decays are also used to determine 
their ratio of branching fractions, which can be used to quantify non-factorizable effects in 
such heavy-to-light decays [9]. 

The paper is organised as follows: the detector is described in Sec. |2j followed by 
the event selection and the relative selection efficiencies in Sec. HI The fit to the mass 
distributions and the determination of the signal yields are discussed in Sec. [4] The 
systematic uncertainties are presented in Sec. [5j and the final results are given in Sec. [6j 



2 Detector and software 

The LHCb detector (TlJ is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r\ < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream. Data are taken with both magnet 
polarities. The combined tracking system has momentum resolution Ap/p that varies 
from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, and impact parameter^] resolution of 20 |J.m 
for tracks with high transverse momentum. Charged hadrons are identified using two 
ring-imaging Cherenkov detectors. 



The trigger 12 consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage which applies a full event reconstruction. 
The events used in this analysis are selected at the hardware stage by requiring a cluster in 
the calorimeters with transverse energy larger than 3.6 GeV. The software stage requires 
a two-, three- or four-track secondary vertex with a high sum of the of the tracks and 
a significant displacement from the primary pp interaction vertices (PVs). At least one 
track should have pp greater than 1.7GeV/c, track fit \ 2 over the number of degrees of 
freedom less than two, and IP \ 2 with respect to the associated primary interaction greater 
than sixteen. The IP \ 2 is defined as the difference between the x 2 from the vertex fit of 
the associated PV reconstructed with and without the considered track. A multivariate 
algorithm is used for the identification of the secondary vertices consistent with the decay 
of a b hadron. 

In the simulation, pp collisions are generated using Pythia 6.4 113] with a specific 



LHCb configuration 14 . Decays of hadronic particles are described by EvtGen 15 



whilst final state radiation is generated using Photos 16 . The interaction of the generated 
particles with the detector and its response are implemented using the Geant4 toolkit 117] 
as described in Ref. |18|. 

2 Impact parameter (IP) is defined as the transverse distance of closest approach between the track and 
a primary interaction. 
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3 Event selection 



The three decay modes, B° — > D~n + , B° — > D~K + and B® — > D~n + , are topologically 
very similar and can therefore be selected using the same event selection criteria, thus 
minimizing efficiency differences between the modes. The B9 % candidates are reconstructed 
from a D7^ candidate and an additional pion or kaon (the "bachelor" particle), with the 
meson decaying to K + tt~it~~ (K + K~tt~). 

After the trigger selection, a loose preselection is made using the B^ a) and D7 s) masses, 
lifetimes and vertex qualities. A boosted decision tree (BDT) (T9J is used to further 
separate signal from background. The BDT is trained on half the B° s — > D~n + data 
sample. The most discriminating variables are the B9^ impact parameter y 2 , the pointing 
angle of the B9 candidate to the primary vertex and the px of the tracks. A cut value for 
the BDT output variable was chosen to optimally reduce the number of combinatorial 
background events, retaining approximately 84% of the signal events. 

The D7, candidates are identified by requiring the invariant mass under the K + tt~tt~ 
(K+K-it-) hypothesis to fall within the selection window 1844 - 1890 (1944 - 1990) MeV/c 2 . 
The relative efficiency of the selection procedure is evaluated for all decay modes using 



simulated events, generated with the appropriate Dalitz plot structures 20 21 . Since the 



analysis is only sensitive to relative efficiencies, the impact of any discrepancy between 
data and simulation is small. 

The final B® — > D~tt + and B° — > D~tt + event samples are obtained after a particle 
identification (PID) criterion, based on the difference in log-likelihood between the kaon 
and pion hypotheses (DLL). A cut on the bachelor particle is placed at DLL(K — 7r)< 0. 
The B° — > D~K + sample is selected by requiring DLL(.K" — tt)> 5 for the bachelor particle. 
The D~ — > K + K~n~ decay is distinguished from D~ — > K + 7t~tt~ decays by imposing 
DLL(K — 7r)> 5 on the kaon candidate with the same charge as the D meson, whilst 
the DLL criteria for the 7r~ and K + are identical between D~ and D~ and are used to 
discriminate D7^ decays from background. The PID performance as a function of and r\ 
of the track is estimated from data using a calibration sample of approximately 27 million 
D*~ — > D (K + 7t~)tt~ decays, which are selected using kinematic criteria only. 



4 Event yields 

The relative yields of the three decay modes are determined from unbinned extended 
maximum likelihood fits to the mass distributions of the reconstructed B9s candidates 
as shown in Fig. [T] In order to achieve the highest sensitivity, the sample is separated 
according to the two magnet polarities, allowing for possible differences in PID performance 
and in running conditions. A simultaneous fit to the two magnet polarities is performed 
for each decay mode, with the peak position and width of each signal shared between the 
two. 

The signal mass shape is described by a Gaussian distribution with power-law tails on 
either side to model the radiative tail and non-Gaussian detector effects. It consists of a 
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Figure 1: Invariant mass distributions of (a) B° — > D~n + (b) B° — > D~K + and (c) 
Bg—t Djii + candidates. The solid line is the result of the fit and the dotted line indicates 
the signal. The stacked background shapes follow the same top-to-bottom order in the 
legend and the plot. The B° s and A b backgrounds in the B° — > D~n + mass distribution are 
invisibly small. The resulting signal yields are listed in Table [TJ For illustration purposes 
the figures include events from both magnet polarities, although they are fitted separately 
as described in the text. 
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and a second, similar but mirrored, function to describe the right tail, resulting in the 
signal mass shape /2Cb(^) = /left(^) + /right (m)- The parameters of the tails are obtained 
from simulated events. The mean \x and the width a of the Gaussian distribution are 
equal in both Crystal Ball functions, and are allowed to vary in the fit. The parameter N 
is a normalisation factor. 

Three classes of background are considered in the fit: fully reconstructed decays where at 
least one track is misidentified, partially reconstructed decays with or without misidentified 
tracks and combinatorial background. The shapes of the invariant mass distributions for 
the partially reconstructed decays are taken from large samples of simulated events. The 
main sources are B° — > D~p + and B° — > D*~7T + (K + ) for the B° — > D~tt + (K + ) sample, 
and 5° ->■ D~ p + and 5° -> L>*~7r+for the 5°-^ D~ir + sample. 

The invariant mass distributions of the misidentified decays are affected by the PID 
criteria. The shapes are obtained from simulated events, with the appropriate mass 
hypothesis applied. The distribution is then reweighted in a data-driven way, according 
to the particle identification cut efficiency obtained from the calibration sample, which is 
strongly dependent on the momentum of the particle. 

Despite the small n — > K misidentification probability of 2.8%, the largest misidentified 
background in the B° — > D~K + sample originates from Cabibbo-favoured B° — > D~tt + 
decays where the bachelor pion is misidentified as a kaon. The shape of this particular 
misidentified decay is determined from data using a high purity sample of B° — > D~tt + 
decays (see Fig. [ija)), obtained by selecting events in a narrow mass window 5200- 
5340 MeV/c 2 . The yield of this prominent peaking background is allowed to vary in the 
fit and is found to be consistent with the expected yield based on the B° — > D~tt + signal 
yield and the misidentification probability. The contamination of B° — > D~tt + events 
in the B® — > D~n + sample can be caused by the misidentification of either pion from 
the D~ decay. The misidentification probability is 2.0% (3.2%) for the higher (lower) 
Pt pion. After selecting the D~ candidate within the mass window around the nominal 
D~ mass |2], the number of misidentified pions is reduced to 0.75% (0.02%). The yield 
of this background is constrained in the fit, based on the B° — > D~tt + signal yield, the 
misidentification probability and their associated uncertainties. 

The yield of A b — > A~TT + dec&ys is allowed to vary in the fit. The cross-feeds from 
B° — > D~K + and B® — > D~n + events in the B° — > D~7t + signal is small, and are 
constrained to their respective predicted yields. In addition, a contribution from the rare 
B° — > D~n + decay is expected with a yield of 3.3% compared to the B®— > D~n + signal, 
and is accounted for accordingly. 
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Table 1: Yields obtained from the fits to the invariant mass distributions. 

Signal Yield 
B ^ D-7T+ 106 197 ± 344 
B°^D-K+ 7664± 99 
B° s ^D;ir + 17419 ±155 



The combinatorial background consists of events with random pions and kaons, forming 
a fake D~~ or D~ candidate, as well as real D~~ or D~ mesons, that combine with a random 
pion or kaon. The combinatorial background is modelled with an exponential shape. 

The results of the fits are presented in Fig. [TJ and the corresponding signal yields are 
listed in Table [TJ The total yields of the decays B° — > D~ir + and B° — > D~K + are used 
to determine the ratio of their branching fractions, while the event yields of the decays 
B° s — > D~7i + and B° — > D~K + are used to measure the average ratio of fragmentation 
fractions. 

The dependence of the relative 6-hadron production fractions as a function of the 
transverse momentum and pseudorapidity of the B9 meson is studied in the ranges 
2.0 < 7] < 5.0 and 1.5 < p T < 40 GeV/c, using D~7r + and 5°-^ D~7r + decays. The 

event sample is subdivided in 20 bins in px and 10 bins in i], with the bin sizes chosen to 
obtain approximately equal number of events per bin. The fitting model for each bin is 
the same as that for the integrated samples, apart from the treatment of the exponent of 
the combinatorial background distribution, which is fixed to the value obtained from the 
fits to the integrated sample. 

5 Systematic uncertainties 

The systematic uncertainties on the measurement of the relative event yields of the 
B° — > D~7i + , B° — > D~K + and B° s — > Djtt + decay modes are related to trigger and offline 
selection efficiency corrections, particle identification calibration and the fit model. 

The response to charged pions and kaons of the hadronic calorimeter used at the 
hardware trigger level has been investigated. As the hardware trigger mostly triggers on 
the high-px bachelor, a systematic uncertainty of 2% is assigned to the ratio of trigger 
efficiencies for the decays B° — > D~K + and B° — > D~tt + , estimated from dedicated studies 
with D*~ — > D°(K + 'k~)'k~ decays. This uncertainty is assumed to be uncorrelated between 
the individual bins in the binned analysis. 

The relative selection efficiencies from simulation are studied by varying the BDT 
criterion, changing the signal yields by about ±25%. The variation of the relative efficiency 
is 1.0% which is assigned as systematic uncertainty. 

The uncertainty on the PID efficiencies is estimated by comparing, in simulated 
events, the results obtained using the D*~ calibration sample to the true simulated PID 
performance on the signal decays. The corresponding uncertainty ranges from 1.0% to 
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1.5% for the different measurements. 

The exponent of the combinatorial background distribution is allowed to vary in the fits 
to the B° — » D~n + and B Q S — > D~tt + mass distributions. By studying D~tt~ and D~K~ 
combinations, it is suggested that the value of the exponent is smaller for the B°— > D~K + 
decays than for the B° — > D~tt + decays, and therefore in the fit to the B° — > D~K + 
candidates the exponent is fixed to half the value found in the fit to the B° — > D~tt + 
sample. The uncertainty on the signal yields due to the shape of the combinatorial 
background is estimated by reducing the exponent to half its value in the fits to the 
B° — > D~7i + and B® — > D~ir + mass distributions, and by taking a flat background for the 
fit to the B° — > D~K + mass distribution. An uncertainty of 1.0% (0.7%) is assigned to 
the relative > D~K + and B° s ^ Dj>k + (B ^ D~7r+) yields. 

The tails of the signal distributions are fixed from simulation due to the presence of 
large amounts of partially reconstructed decays in the lower sidebands. The uncertainty on 
the signal yield is estimated by varying the parameters that describe the tails by 10%. The 
uncertainty from the shape of the central peak is taken from a fit allowing for two different 
widths for the Crystal Ball functions in Eq. [2| leading to a 1.0% (0.8%) uncertainty on the 
relative > D~K + and B° s ^ D;n + (B ^ D~<k + ) yields. 

The contribution of charmless B decays without an intermediate D meson is ignored 
in the fit. To evaluate the systematic uncertainty due to these decays, the B mass spectra 
for candidates in the sidebands of the D mass distribution are examined. A contribution 

Table 2: Systematic uncertainties for the measurement of the corrected ratio of event yields 
used for the measurements of f s /fd and the relative branching fraction of B°— > D~K + . 
The systematic uncertainty in and t] bins is shown as a range in the last column, and 
the total systematic uncertainty is the quadratic sum of the uncorrelated uncertainties. 
The systematic uncertainties on the ratio of B° — > D~n + and B° s ^- D~n + yields that are 
correlated among the bins do not affect the dependence on or rj, and are not accounted 
for in the total systematic uncertainty. 



Source i£g=& (%) (%) 



Detector acceptance 








and reconstruction 


0.7 


0.7 


2.0-2.9 


Hardware trigger efficiency 


2.0 


2.0 


0.8 


Offline selection 


1.2 


1.1 


1.2 


BDT cut 


1.0 


1.0 


1.5 


PID selection 


1.0 


1.5 


1.1 


Comb, background 


0.7 


1.0 


0.8 


Signal shape (tails) 


0.5 


0.6 


[correl.] 


Signal shape (core) 


0.8 


1.0 


[correl.] 


Charmless background 


0.4 




[correl.] 



Total 3.1 3.4 3.2-3.8 
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of 0.4% relative to the signal yield is found in the B° — > D~tt + decay mode, and no 

contribution is seen in the other modes. For the B° — > D~n + decay mode no correction is 

applied and the full size is taken as an uncertainty. No systematic uncertainty is assigned 

for the other decay modes. 

The various sources of the systematic uncertainty that contribute to the uncertainties 

on the ratios of signal yields are listed in Table |2| No uncertainty is associated to the 
o 1 — 1 

A b — > /l^vr+background, as the yield is allowed to vary in the fit. Other cross checks, 
like varying the B° — > D~n + yield in the B° s — > D~tt + fit or including A b — > A~K + in the 
B° — y D~K + fit, show a negligible effect on the signal yields. 

All systematic variations are also performed in bins, and the corresponding relative 
changes in the ratio of yields have been quantified. Variations showing correlated behaviour 
do not affect the slope and are therefore not considered further. 



6 Results 

The relative signal yields of the decays B° — > D~n + , B° — > D~K + and B® — > D~n + are 
used to determine the branching fraction of the decay B° — > D~K + , and the ratio of 
fragmentation fractions f s /fd ■ 

The efficiency corrected ratio of B° — > D~K + and B° — > D~tt + signal yields results in 
the ratio of branching fractions 

B(B°^D-K+) . . , 

— \— ^_ J = 0.0822 ± 0.0011 (stat) ± 0.0025 (syst). 

B [B u — > D 7T+) 

This is combined with the world average branching fraction B (B°— > D~n + ) = (26.8 ± 
1.3) x 10~ 4 |2], to give 

B (B°^ D-K + ) = (2.20 ±0.03 ±0.07 ±0.11) x 10" 4 , 

where the first uncertainty is statistical, the second is systematic and the last is due to 
the uncertainty on the B° — > D~n + branching fraction. 

The ratio of fragmentation fractions is determined from the efficiency corrected event 
yields. The ratio of efficiencies is 0.913 ± 0.027. This results in 

f 1 
(0.261 ± 0.004 ± 0.017) x 



fd MoNf 

= 0.238 ±0.004 ±0.015 ±0.021, 

where the first uncertainty is statistical, the second is systematic containing the sources 
listed in Table [2] as well as errors from external measurements, and the third is theoretical, 
due to the knowledge of Af a and Afp. The last source is dominated by the uncertainty on 
the form factor ratio. 

This measurement supersedes and is in agreement with the previous determination with 
hadronic decays [3]. It also agrees with the previous measurement based on semileptonic 
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Figure 2: Ratio of fragmentation fractions f s / fa as functions of (a) pr and (b) 77. The 
errors on the data points are the statistical and uncorrelated systematic uncertainties 
added in quadrature. The solid line is the result of a linear fit, and the dashed line 
corresponds to the fit for the no-dependence hypothesis. The average value of pp or 77 
is determined for each bin and used as the center of the bin. The horizontal error bars 
indicate the bin size. Note that the scale is zero suppressed. 



decays El. The two independent results are combined taking into account the various 
sources of correlated systematic uncertainties, notably the D7, branching fractions and 
B9 s lifetimes, to give 

f 

^ = 0.256 ±0.020, (3) 

Id 

which supersedes the previous measurement from LHCb. 

The value of f s /fd in bins of j?t or 77 is determined using the B° s — > D~tt + and 
B° — > D~tt + decay modes and is presented in Fig. [2] A linear \ 2 fit gives 



fs/fd M = (0.256 ± 0.020) + (-2.0 ± 0.6) x IO- 3 / GeV/c x (p T - (p T » 
fs/fd (V) = (0-256 ±0.020) + (0.005 ± 0.006) x (77- (77)), 

with (px) = 10.4 GeV/c and (77) = 3.28. The data points are normalised with a scale 
factor to match the average value of 0.256. The uncertainty associated to this parameter 
is taken from Eq. [3j whilst the error from the fit is 0.003 for both and r). 

The p- value for this linear fit is found to be 0.16 (0.87) for (if). The observed slope 
for the dependence on the transverse momentum of the B? s n meson deviates from zero 
with a significance of three standard deviations. No indication of a dependence on rj(B) is 
found. 
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7 Conclusions 



The relative production rate of B® and B mesons is determined using the hadronic decays 
B° s ^ DJ7T+ and B° — > D~K + resulting in f s /f d = 0.238 ± 0.004(stat) ± 0.015(syst) ± 
0.021(theo). This value is consistent with a previous LHCb measurement based on 
semileptonic decays, with which it is averaged to obtain f s /fd = 0.256 ± 0.020. The ratio 
of fragmentation fractions f s /fd is determined as a function of the transverse momentum 
and pseudorapidity of the B®^ meson, and a variation consistent with a linear dependence 
on the transverse momentum of the the B9\ meson is observed, with a significance of 
three standard deviations. In addition, the ratio of branching fractions of the decays 
B ^ D~K+ and B° D--K+ is measured to be 0.0822 ± 0.0011 (stat) ± 0.0025 (syst). 



Acknowledgements 

We express our gratitude to our colleagues in the CERN accelerator departments for 
the excellent performance of the LHC. We thank the technical and administrative staff 
at the LHCb institutes. We acknowledge support from CERN and from the national 
agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); NSFC (China); CNRS/IN2P3 
and Region Auvergne (France); BMBF, DFG, HGF and MPG (Germany); SFI (Ireland); 
INFN (Italy); FOM and NWO (The Netherlands); SCSR (Poland); ANCS/IFA (Romania); 
MinES, Rosatom, RFBR and NRC "Kurchatov Institute" (Russia); MinECo, XuntaGal 
and GENCAT (Spain); SNSF and SER (Switzerland); NAS Ukraine (Ukraine); STFC 
(United Kingdom); NSF (USA). We also acknowledge the support received from the ERC 
under FP7. The Tierl computing centres are supported by IN2P3 (France), KIT and 
BMBF (Germany), INFN (Italy), NWO and SURF (The Netherlands), PIC (Spain), 
GridPP (United Kingdom). We are thankful for the computing resources put at our 
disposal by Yandex LLC (Russia), as well as to the communities behind the multiple open 
source software packages that we depend on. 



References 

[1] LHCb collaboration, R. Aaij et al., First evidence for the decay B® — > Phys. 



Rev. Lett. 110 (2012) 021801, [arXiv: 1211 . 2 674 



[2] Particle Data Group, J. Beringer et al, Review of particle physics, Phys. Rev. D86 
|(2012) 010001) 

[3] LHCb collaboration, R. Aaij et al., Determination of f s /fd for 7 TeV pp collisions 
and measurement of the B° — > D~K + branching fraction, Phys. Rev. Lett. 107 (2011) 
|211801| |arXiv : 1 106 . 4435| 

[4] LHCb collaboration, R. Aaij et al., Measurement of b hadron production fractions in 
7 TeV pp collisions, |Phys. Rev. D85 (2012) 032008| |arXiv: 1111 .2357[ 



10 



[5] R. Fleischer, N. Serra, and N. Tuning, New strategy for B s branching ratio measure- 



ments and the search for new physics in B° s — > , Phys. Rev. D82 (2010) 034038 
larXiv: 1004 .39821 

[6] Heavy Flavor Averaging Group, Y. Amhis et ai, Averages of b-hadron, c-hadron, and 
r-lepton properties as of early 2012, arXiv7l207 . 1 158} 



[7] CLEO collaboration, S. Dobbs et ai, Measurement of absolute hadronic branching 
fractions of D mesons and e + e~ — >■ DD cross sections at the ?/>(3770), |Phys. Rev] 
|D76 (2007) 11200lHarXiv: 0709 3783] 



CLEO collaboration, J. Alexander et ai, Absolute measurement of hadronic branching 
fractions of the D+ meson, |Phys. Rev. Lett. 100 (2008) 161804] |arXiv:0801 .0680[ 

[9] R. Fleischer, N. Serra, and N. Tuning, Tests of factorization and SU(3) relations in B 
decays into heavy-light final states, Phys. Rev. D83 (2011) 014017[ |arXiv: 1012 . 2784| 

[10] J. A. Bailey et ai, B s — > D s /B — > D semileptonic form-factor ratios and their 
application to BR(B° S ->• [Phys. Rev. D85 (2012) 114502HarXiv: 1202.6346} 

Erratum-ibid. D86 (2012) 039904. 



[11] LHCb collaboration, A. A. Alves Jr. et ai, The LHCb detector at the LHC, JINST 3 
|(2008) S08005| 



[12] R. Aaij et al, The LHCb trigger and its performance, arXiv: 1211.3055 



[13] T. Sjostrand, S. Mrenna, and P. Skands, PYTHIA 6.4 physics and manual, JHEP 05 
|(2006) 0~26l|arXiv:hep-ph/0603175 



[14] I. Belyaev et al., Handling of the generation of primary events in Gauss, the LHCb 



simulation framework, Nuclear Science Symposium Conference Record (NSS/MIC) 
|IEEE (2010) 1155 



[15] D. J. Lange, The EvtGen particle decay simulation package, Nucl. Instrum. Meth 
|A462 (200iyi52j 

[16] P. Golonka and Z. Was, PHOTOS Monte Carlo: a precision tool for QED corrections 
Z and W decays, Eur. Phys. J. C45 (2006) 97\ |arXiv : hep-ph/0506026] 



in 



[17] GEANT4 collaboration, J. Allison et ai, Geant4 developments and applications, 
|IEEE Trans. Nucl. Sci. 53 (2006) 270[ GEANT4 collaboration, S. Agos tinelli et al, 
GEANT4: A simulation toolkit, |Nucl Instrum. Meth. A506 (2003) 250| 

[18] M. Clemencic et ai, The LHCb simulation application, Gauss. : design, evolution and 
experience, |J. of Phys. : Conf. Ser. 331 (2011) 032023| 

[19] L. Breiman, J. H. Friedman, R. A. Olshen, and C. J. Stone, Classification and 
regression trees, Wadsworth international group, Belmont, California, USA, 1984. 



11 



[20] BaBar collaboration, P. del Amo Sanchez et ai, Dalitz plot analysis of — > 
K+K-1T+, | |Phys. Rev. D83 (2011) 0520011 |arXiv: 1011 .4190| 

[21] CLEO collaboration, G. Bonvicini et ai, Dalitz plot analysis of the D + — > K~tt + tt + 
decay, Phys. Rev. D78 (2008) 052001[ |arXiv:0802.4214[ 

[22] T. Skwarnicki, A study of the radiative cascade transitions between the Upsilon-prime 
and Upsilon resonances, PhD thesis, Institute of Nuclear Physics, Krakow, 1986, 
IDESY-F31-86-021 



12 



